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pLn ; Abstract 

In this paper we carry out a complete classification of the Lie point symmetry groups associated 
with the quadratic Lienard type equation, x + f{x)x 2 + g(x) = 0, where f{x) and g{x) are arbitrary 
functions of x. The symmetry analysis gets divided into two cases, (i) the maximal (eight parame- 
ter) symmetry group and (ii) non-maximal (three, two and one parameter) symmetry groups. We 
identify the most general form of the quadratic Lienard equation in each of these cases. In the 
case of eight parameter symmetry group, the identified general equation becomes linearizable as 

m ■ 

(f-) ' well as isochronic. We present specific examples of physical interest. For the nonmaximal cases, 

o ' 

the identified equations are all integrable and include several physically interesting examples such 
as the Mathews-Lakshmanan oscillator, particle on a rotating parabolic well, etc. We also analyse 

j— i ' 

the underlying equivalence transformations. 
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I. INTRODUCTION 



Ordinary differential equations (ODEs), especially nonlinear ones, are very useful in the 
formulation of fundamental natural laws and technological problems for a long time. The 
Lie algebraic properties of these equations are one of the basic aspects which got attention 



after Lie's initial work in which 
can be obtained from his theory 



re discovered that all the integration methods for ODEs 



-llOj. Lie gave a classification of ODEs in terms of their 
symmetry groups, thereby identified the full set of equations which could be integrated or 
reduced to lower order equations by his method. Since then several contributions have been 
made on symmetry group classification of ODEs 
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17(| . In particular it has been shown 



that any second order nonlinear ODE which admits eight parameter Lie point symmetries 
is linearizable to free particle equation through point transformations lfl. 

Recently, Pandey et al. [12j have studied the Lie point symmetry properties of a general 
Lienard type equation 

x + f(x)x + g(x) = 0, (- = ^) (1) 



where f(x) and g(x) are arbitrary smooth functions of x and overdot denotes differentiation 
with respect to t, and identified several interesting integrable and linearizable equations. 
They divided their analysis into two parts. In the first part they isolated equations that 
admit lesser parameter Lie point symmetries, and in the second part identified equations 
that admit maximal (eight parameter) Lie point symmetries. They proved the integrability 
of all the equations obtained in the first part either by providing the general solution or 
by constructing a time independent Hamiltonian. In the second part they discussed the 
linearizing transformations and solutions for all the nonlinear equations identified under 
this category. 

Yet another general second order nonlinear differential equation which is of high physical 



and mathematical interest 



id, 
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19] is of the form 



A(x, x, x) = x + f(x)x + g(x) = 



(2) 



where f(x) and g(x) are arbitrary functions of x and overdot denotes differentiation with 
respect to t. We will designate Eq. (j2j) as a quadratic Lienard type system for convenience 
(corresponding to x 2 term in (j2J)). For example, the one dimensional Mathews-Lakshmanan 

2 

(ML) oscillator with f(x) = — 14 x f 2 and g(x) = ,T 2 belongs to this class 22]. It exhibits 
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only one Lie point symmetry even though it has been proved to be linearizable with the 



help of nonlocal transformations. Due to its unusua' 
authors at the classical as well as quantum levels [23 



property it has been studied by many 
24J | . Another example is the motion of 

27 L Several 



a particle on a rotating parabolic well with f(x) = rrf 2 and g(x 



i+\x 2 aw 1+Ax 2 
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theorems on the isochronous cases belonging to (J2J) exist in the literature 18l42l|. Even 
though the integrability properties of the general form of Eq. (j2J) has been discussed in the 
literature and it can be integrated in the form of a quadrature, its general properties have 
been hardly studied. In this paper, we systematically identify and classify all the equations, 
belonging to class (J2J), which admit one, two, three and eight parameter symmetry groups, 
from a group theoretic point of view and explore certain interesting properties associated 
with them. 

The main objective of this paper is to carry out a detailed Lie point symmetry analysis of 
Eq. (j2J). In this way we study the linearizable, and integrable (but not linearizable by point 
transformations) cases separately. Firstly, we consider the linearizable case and find the 
general form of Eq. (j2J) for which it admits eight point symmetry generators. The general 
form of Eq. (j2J) in this case turns out to be 

where g\ and #2 are constant parameters. Apart from this we also show that the system (j2J) 
additionally follows the isochronous condition, 

9x + f9 = 9u (4) 

where g\ is a constant and subscript denotes differentiation with respect to x. We prove the 
significant result that (j4]) implies (jSJ). We also discuss some specific examples of physical 
interest belonging to this class. 

Secondly, we consider the integrable cases of Eq. (j2J) with lesser parameter symmetries. 
The general form of the equations which show three and two parameter symmetry generators, 
respectively, are as follows: 

(i) x + f(x) x 2 + gi e- 1 fdx (Ai + J fdx dxj ' = 0, (5) 

x + f{x)x 2 + g 1 e-f fdx (\ 1 + J ef fdx dx^ " = 0, A 2 ± 1, 4, (6) 
where g±, X± and A 2 are constants. 



The rest of the cases of Eq. (jSj) correspond to one parameter Lie point symmetry group. 

2 

In fact, ML oscillator with f(x) = — 1H A A a: 3 ,2 > d( x ) — i+az 2 belongs to this case. To ensure the 
existence of integrability of such systems, one can associate more general symmetries such 
as the A-symmetries 28l. l29|. 

The plan of the paper is as follows. In Sec. II we deduce the determining equations for 
the infinitesimal symmetries. The general form of the equation is obtained for the maximal 
symmetry group in Sec. Ill and its symmetry generators corresponding to sl(3, R) algebra 
are deduced. The isochrnocity of the system is verified in Sec IV. Some physically interesting 
examples of this class of equations are discussed in Sec V. Sec. VI deals with the nonmaximal 
case and the associated equations are derived here. The symmetry property with f(x) or 
g(x) as zero is discussed in Sec. VII. The equivalence transformations are investigated in 
Sec. VIII. Finally, conclusion is given in Sec. IX. 



II. DETERMINING EQUATIONS FOR THE INFINITESIMAL SYMMETRIES 

We consider the quadratic Lienard type system of the form (J2J). We will assume f(x) ^ 
and g(x) 7^ to start with. Let Eq. (j2J) be invariant under an one parameter group of 
symmetry transformations, 

T = t + e£(t,x) + 0(e 2 ) 

X = x + er]{t,x) + 0(e 2 ), e < 1. (7) 



An operator G given by 



c) 

G(t,x)=£(t,x) — + 7i(t } x)— (8) 



is said to be an infinitesimal generator of the one parameter Lie point symmetry group of 
transformations for Eq. (j2J) iff 

G^\A) U =0 =0, (9) 



or equivalently 



where 



BA dA m dA (2] dA> 

£ « + '& + ' m » + ' m «)i*-= " (10) 



G « = G (l) + I) (2)iL G <1> = G + „(1)JL 



dt dt 1 dt dt v ' 

and 

d d , d . . 

Substituting Eqs. ffTTj) and f|T2|) in Eq. ffTUj) and equating different powers of x m , m = 
0, 1, 2, 3, to zero, we obtain the determining equations, 

U - = 0, (13) 
V xx + fVx+Vfx-2£t x = 0, (14) 
2Vt* + 2fVt -£tt + 3gZx = 0, (15) 

mt + v9x- gvx + 2 g£t = o, (16) 

where subscripts denote partial derivatives. Solving Eq. (TT3"j) we get 

e = 6(t)9f(x) + a(t), (17) 

where 3(x) = J F(x)dx and F(x) = e^^ dx . Here a(t) and 6(t) are arbitrary functions of 
t. Substituting £ tx = bF from (fTTl) into ffT4|) and then integrating it twice with respect to x, 
we get 

7] = 2bG 1 (x) + c(t)G 2 (x) + d(t)G 3 (x), (18) 

where Gi(x) = - ^M^)^ ; G 2 (x) = G 3 (x) = -^^y and c(i) and d(t) are arbitrary 

functions of t. Now with the above forms of £ and 77, Eqs. ( 1T51) and ( fl6l) can be written as 

2/(2&Gi + cG 2 + dG 3 ) +4 6Gi x + 2cG 2 x + 2dG 3iB -63-S + 3#6F = 0, (19) 

and 

g x (2 6 Gi + c G 2 + d G 3 ) - g (2 6 G lx + c G 2x + d G 3x - 2 d - 2 6 Of) 

+26'Gi + cG 2 + dG 3 = 0. (20) 

In analyzing the above system of Eqs. (fl9~|) and (120"]) . we can distinguish two separate 
cases. 

(i) Case 1 : b ^ : In this case we can write Eq. ( fl9l) as 

* = sTf • (21) 
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Substituting /, the above form of g and Gi s,i = 1,2,3, into (120]) and equating the coef- 
ficients of various independent functions of x, a set of determining equations can be ob- 
tained. Solving the resultant determining equations one can get the symmetry functions 
a(t), b(t), c(t) and d(t) which in turn fix the form of the associated symmetries. 
(ii) Case 2 : 6 = : In this case Eqs. (IT51) and f[2"0"j) are simplified to 

2f(cG 2 + dG 3 ) + 2 (cG 2x + dG 3x ) - a = 0, (22) 

and 

g x (cG 2 + dG 3 )-g (cG 2x + dG 3x - 2d) + cG 2 + dG 3 = 0. (23) 

Substituting the forms of G 2 and G 3 in Eqs. (|22|) and (1231) one can get the form of g as well as 
the determining equations for the symmetry functions. Solving the associated determining 
equations, one can identify the corresponding symmetries. 

Consequently, we will investigate the two cases, (i) b ^ and (ii) 6 = 0, separately in 
sections III-V and VI, respectively, and show that the b ^ case corresponds to maximal 
(eight parameter) symmetry group, while the 6 = case corresponds to nonmaximal (three, 
two, one parameter) symmetry group of transformations. 

III. GENERAL FORM OF THE EQUATION FOR 6/0 CASE - EIGHT PARAM- 
ETER SYMMETRIES 

We now consider the case 6 ^ in Eqs. (|T9j) and (120]) . For the sake of generality we 
consider / = f(x). Now substituting the values of Gi s,i = 1,2,3, (given below Eq. (|T8l) ) 
and their derivatives in terms of f(x) in Eq. (1211) we arrive at 

g = g ie -f f(x)dx J ef f{x)dx dx + g 2 e-f f{x)dx , (24) 

where g± — — | = constant and g 2 = — = constant as g is a function of x only and is 
free from t. 

Now for the above form of g(x) Eq. (j2J) can be written as 

x + f{x) x 2 + g x e~ 1 f(x)dx J e f f ^ dx dx + g 2 e~ * f{x)dx = 0. (25) 

To explore the symmetry group of the above Eq. (|25|) . we substitute the values of Gi s 
and their derivatives as well as g in (120|) . and equate the resultant coefficients of various 
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independent functions of x to zero. Consequently we get the following set of determining 
equations for the functions ait), b(t), c(t) and d(t), that is 

b + gib = 0, d — 2c — 3g 2 b = 0, 
c + 2 gi a = 0, d + g\ d + 2 g 2 a — g 2 c = 0. (26) 

Note that g\ and #2 are system parameters fixed by the form of g[x) in ff2^|) . Solving the 
above system of equations consistently, we can express the form of the functions a(t), b(t), c(t) 
and d(t) as 

a (t) =ai- — (c 2 sin (2 v /^t) + c 3 cos (2 v /^"t) -2g 2 h sin (y^t) -2g 2 b 2 cos (^t)) , 
6(t) = 6i sin (v^It) + & 2 cos (y/glt) , 

c (t) = ci - 1 (c 2 cos (2 y/glt) - c 3 sin (2 y/glt) -Ag 2 b 1 cos (x/^l*) + 4 # 2 &2 sin {y/g[t)) , 

(2 (t) = di sin (v^T*) + d 2 cos (v^i*) H jm ( 2 ciy^i ~ c 2 cos (2 y^t) + c 3 sin (2 y^)) , 

(27) 

where ai, &i, 62, ci, C2, c 3 , a?i and are arbitrary constants, which indeed constitute the 
eight symmetry parameters corresponding to the maximal symmetry group. Thus we can 
conclude that the specific differential Eq. (1251) admits the maximal symmetry group. The 
corresponding infinitesimal symmetries (which follow from ( TT71) and ( TT8l) ) are 

sin (y/g~it) (gi J e^^ dx dx + g 2 ) cos (yfg~[t) {g x J e^^ dx dx + g 2 ) 



i = ai + bi 



9i 9i 



77 = 61 



sin (a/2 git) cos (a/2 git) 

e -ff(x)dx (J e ff(x)dx dx ^ CQS (^t) ( 5l J e ff(*)d*dx + 2 2 ) 



-6 



'01 

e -ff(x)dx y e ff(x)dx dx ^ gin (^t) (^ J e tt^ dx dx + 2 g 2 ) 



2 — 

e -ff(x)dx( gi j e ff(x)dx dx + ^ e -ff(x)dx CQS J 2 J e If(x)dx dx + ^ 
+ Cl C2 

5i 2of 2 
+ c 3 e-^sm(2^)( f/ e^)^x + , 2 ) + ^ _ 

+rf2Cos( v ^Tt)e-^^. (28) 
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The associated infinitesimal generators read 



X, 



X-, 



x 4 
x 5 



d_ 

dt' 

sin (y/fht) (9i J ef f{x)dx dx + g 2 ) 8 



9i dt 
e - Jf(x)dx y e jf(x)dx dx ^ CQS (^gfr) (gx j e$f^ dx dx + 2g 2 ) d 



gl dx' 
cos (Jg{t) (gi J effW^dx + g 2 ) d 
9i dt 
e - Jf(x)dx y e ff(x)dx dx j gin ygit) (gx j ef f ^ dx dx + 2g 2 ) d 



Jgi dx' 
e - M*)dx( gi J e ff(x)dx dx + Q 

9i dx' 
sin (y/Tglt) d e~ f f(x)dx cos (2 y/git) {gx J e^^ dx dx + g 2 ) 8 



2g x dt 2gl /2 dx' 

cos (y/Zgit) d e- ffW dx sin (2 y/git) {gx J e$ f ^ dx dx + g 2 ) d 
6 2gx dt 2g^ 2 dx' 

X 7 = S in(y/g-xt)e-f^ dx -^, 

X 8 = sin(^It)e-/^^. (29) 
One can easily check that these eight generators lead to the sl(3, R) algebra. 



IV. ISOCHRONOUS CONDITION, LINE ARIZ ABILITY AND THE NATURE OF 
SOLUTION OF EQ. ([25ft 

A. Isochronicity condition 

A dynamical system is called isochronous if it features in its phase space an open, fully- 
dimensional region where all its solutions are periodic in all its degrees of freedom with the 
same, fixed, period [3(1 ] . The linear harmonic oscillator is the prototype of an isochronous 
system and all other isochronous systems are isoperiodic with the harmonic oscillator. Eq. 
([2]) can be mapped on to the linear harmonic oscillator, 

X + u)%X = 0, (30) 



S 



with the following invertible transformation 

X = h(x), (31) 

where h(x) is a function of x alone, provided for suitable choice of f(x) and g(x) Eq. (j2j) is 
isochronic. Now substituting the derivatives of X in equation fl30|) we get 



X+ h'(x) X +UJ °h'(x)~^ [6Z) 

where prime denotes differentiation with respect to x. Comparing Eq. (|32j) with (|2j) we get 
the following two relations 

<*£g =«(*>■ (34) 



From Eq. ( 133]) one can get 



%) = /n y e$ f{x)dx dx + /iz, (35) 



where /ii and h% are integration constants. Now substituting ( 135]) and its first derivative 
into Eq. ( I34"j) we can obtain 

h(x) = ^ g(x)e^^ dx . (36) 
Comparing the forms of h(x) from Eq. (135]) and (|36l) we obtain the condition 

fti 3(3) f{x)dx = ul (h x J ef f{x)dx dx + h^j . (37) 

Differentiating both sides of ( )37|) with respect to x and integrating the resultant equation 
we find 

g = g l6 - f(x)dx J ef f{x)dx dx + g 2 e-^ f(x)dx , (38) 

where has been replaced by gi for analogy. Eq. ([581 is exactly the same as the one 
we obtained with the help of symmetry method, see Eq. ( 1241) . Now inverting the relation 
X = h(x) in Eq. (131 p . we find that for the above form of g, the solution is isochronous. 
Thus we can conclude that the periodic solutions of Eq. (125]) are all isochronous. 
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B. Linearizability Condition 



It has been proved that the linearization of a scalar second order ODE, x + f(t, x, x) = 0, 
via point transformation has the cubic in first derivative that is 

x = P(t, x)x 3 + Q(t, x)x 2 + R(t, x)x + S(t, x), (39) 

with the coefficients A, B, C and D satisfying the following two invariant conditions, 

3P tt + 3P t R - 3P X S + 3PR t + R xx - 6PS X + QR X - 2QQ t - 2Q tx = 0, 
6P t S - 3Q X S + 3PS t + Q tt - 2R tx - 3QS X + 3S XX + 2RR X - RQ t = 0, (40) 

where the suffix refers to partial derivatives. On comparing §2§ with (I3"9"|) . we find the 
condition of linearizability from equation (1401) as 

9xx + fg x + 9 fx = 0. (41) 

The integration of the above equation gives the form of g(x) as 

g = g^-ffto** J e f f{x)dx dx + g 2 e-f f{x)dx . (42) 

So we see that the linearizability criterion also suggests the same form of g as the one we 
obtained from symmetry method and by isochronicity condition. This is in conformity with 
the fact that the system (|25p admits eight Lie point symmetries as shown above. 

V. SPECIAL CASES OF MAXIMAL SYMMETRY GROUP 

In this section we consider some physically interesting special cases by fixing the form 
of f(x) in Eq. ([25]) . The determining equations will be the same as given by Eq. ( j26|) 
and hence the forms of the functions a(t) , b(t) , c(t) and d(t) will also be the same. Then 
their symmetries can be obtained directly from Eq. (|2"5|) for appropriate choice of f(x) and 
parameters, which leads to the infinitesimal generators. 



A. f(x) = A = constant 

Fixing the function f(x) as a constant (= A), we get the form of Eq. (T5]) from Eq. ff25l) as 
a generalized Morse oscillator, 

x + Xx 2 + ^(l-e- Xx ) = 0. (43) 
A 
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Note that in the limit A — > 0, Eq. (143|) reduces to the linear harmonic oscillator 

x + u^x = 0. (44) 

The form of £ and r\ associated with Eq. (l4*3|) will be obtained directly from equation ( 128]) 
and can be written as 

sin (w t) (e Ax - 1) cos (w t) (e Ax - 1) sin (2 w t) cos (2 w t) 

K (£, X) — di + Ol r t>2 7 C 2 g C 3 7, 2 ' 

A A Z Wg L Wg 



, , wo cos (w t) (e Xx -2) w sin (w t) (e A *- 2) I - e~ Xx 
V(t,x) = h — b 2 — + ci 

cos (2 w t) (1 - e- Xx ) sin (2 ui t) (1 - e~ Ax ) 

_C 2 7T\ 1" c 3 7f\~ 

2Aw 2 A w 

+dx sin (w t) e~ Aa; + d 2 cos (w t) e~ Ax - (45) 

The corresponding symmetry generators thus turn out to be 

d _ sin (w t) (e Xx - 1) d w cos (w t) (e Xx - 2) 5 

1 ~ at ' 2 ~ a at + A2 a? 

cos (w t) {e Xx - 1) 5 w sin (w t) (e Aa; - 2) a 1 - e- Ax a 



-^3 — ; 777 T5 77- ) X. 



A at A 2 ax' " 4 A dx 



sin (2 w t) a _ cos (2 w t) (1 - e~ Xx ) d 

5 2 Wq at 2 A w dx 1 
x _ cos (2 w t) a + sin (2 w t) (1 - e~ Aa; ) a 

6 2u)n dt 2 A w arc ' 



X 7 = sin (w t) e~ Xx —, X 8 = cos (w t) e" Aa; — . (46) 

One can write down the solution of equation ()43p by transforming it to the linear harmonic 
oscillator equation as mentioned above and it is found to be 

x = ~ln(l - AAsin(w t + 6)), < A < -. (47) 
A A 

Here A and 8 are arbitrary constants. For A < 1 we have physically acceptable (isochronous) 
periodic solutions with frequency exactly the same as that of the linear harmonic oscillator 
case, A = 0. The Hamiltonian of equation ()43l) can be then written as 

H=^e-^+ 1 -ul(l-e Xx r, (48) 

where p = X 2 e 2Xx x is the conjugate momentum. 

The phase portrait underlying Eq. ( l4*3l) corresponding to the Hamiltonian (I48p for four 
different sets of A values with A = 0.8, w = 2.6 and 8 = is shown with four different colors 
in Fig.l. Note that regular motion is restricted to —00 < x < jln2. Outside this region, 
the solution becomes singular periodically. 
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FIG. 1: (Color online) Phase portrait of equation (j43]) for different values of A corresponding to 
the Hamiltonian (1481) 



If we consider the form of the function f(x) as / = — 1 2 ^ x , Eq. (1251) becomes 

2 A 

x — x 2 + Unix + \x 2 ) = 0. (49) 

1 + Ax u 

Note that in the limit A — > 0, Eq. (H9j) reduces to that of a linear harmonic oscillator. The 
corresponding infinitesimal symmetries associated with ( 141?]) can be directly written with the 
help of Eq. fl28]) as 

x sin (co>o^) x cos (oj t) sin (2 u; t) cos (2 oj £) 

= ai + bl ^Txx- + h ^Txx-- C2 ^r- C3 ^^^ 



(x + A x 2 ) cos (2 a; £) (x + A x 2 ) sin (2 oo t) 

— C2 T-j 1" C3 ^ 

2 cj 2 o; 

+di sin (u Q t) (1 + A x) 2 + d 2 cos (w *) (1 + A x) 2 . (50) 

The corresponding symmetry generators are given by 

d xsin(u;ot) d | u cos (a; £) (1 + 2 A x) (9 



dt 1 1 + Xx dt A 2 9x ' 



A 

sin (2 wot) 9 (x + Ax 2 ) cos (2u) t) d 



xcos(o> t) d ujq sin (a> £) (1 + 2 A x) <9 2 \ ^ 

* 3 = 1 + Ax dt A2 ^4 = (x + Ax)-. 





2uq dt 2 w <9x' 
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cos (2 u t) d (x + Ax 2 ) sin (2u t) d 
dt 2 dx ' 



2u 2 



X 7 = sm(u t) (1 + Xx) 2 — , X s = cos{u Q t){l + \x) 2 —. 



(51) 



Equation fH9|) can be transformed to the linear harmonic oscillator equation as discussed in 
Sec. IV and the isochronous solution can be then written as 



1 - XAsm(u t + 5) ~ X 



(52) 



where uq and 5 are constants. For A < j-, we have physically acceptable periodic solutions. 
The frequency of oscillation is exactly the same as that of the linear harmonic oscillator 
case, A = 0. The Hamiltonian of fl49l) can be written as 



(53) 



where the conjugate momentum is p 



The phase portrait of Eq. 



for four 



(i+A in- 
different sets of A values with A = 0.8, u = 1.5 and 5 = is shown with four different colors 

in Fig.2. Note that x is restricted to the region — i < x < oo for periodic solutions without 

singularity. 




-0.5 0.0 0.5 1.0 1.5 2.0 



FIG. 2: (Color online) Phase portrait of equation (I49p for different values of A corresponding to 
the Hamiltonian (1531) 
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VI. NON-MAXIMAL SYMMETRY: CASE 6 = 



To explore the forms of the equations having lesser parameter Lie point symmetries we 
consider the case 6 = vide Eqs. ffTU]) and f[2"Uj) . Eq. (1221) provides condition between 
the symmetry functions a(t), c(t) and d(t), whereas Eq. (j2"51) determines the form of g 
corresponding to the given form of /. 

Let us now define a function M = 2 (cG 2 + dG$). Then Eq. (I22l) can be rewritten as 

/M + M x - a = 0. (54) 

Integrating Eq. (J54"j) we get 

MF = a$f(x) + d 1 , (55) 
where d± is an integration constant. Consequently, we have 

fx, 7 

2cG 2 + 2dG 3 = a- + -i (56) 

that is 

(2c-a)G 2 + (2d-di)G 3 = 0. (57) 

As G 2 7^ and G 3 ^ 0, which are two distinct functions of f(x), we are left with two 
constraints 

2 c = a, 2d = d 1 . (58) 



Next we note that Eq. ( 1231) is a first order differential equation in g. For a given form 
of / in (|23|) one can get g which will then decide the form of (j5J- Rewriting Eq. ( )23|) with 
d = (from (158|) ). we have 

9 *- 9 \ cG 2 + dG 3 ) + cG 2 + dG 3 =0 - (59) 
Since g should be a function of x alone (vide Eq. (j2J)), we choose 

d . 2d c 

- = Ai, — = A 2 , -- = A 3 , 60 

c c c 

where Ai, A 2 and A 3 are constants. Then Eq. (1591) can be rewritten as 

/ Got ~\~ Ai G 3t — Ao \ A 3 (jro , , 

9*-9{ c ' r ) + r l xr =0- (61) 

\ Lr 2 -+- AiLr 3 / Lt 2 -+- Ai<_t 3 
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Analyzing the relations 2 c = a (from (EE}) and A 2 = — (from (l6"0| )). we get 

(A 2 - 4)c = 0, (62) 

provided (G 2 + A1G3) 7^ or CG2 + (IG3 7^ 0. Hence we are left with three different 
possibilities as (i) X 2 = 4, c ^ and [cG^ + (IG3) 7^ 0, (ii) A2 7^ 1, 4, c = and (cG*2 + 
dGs) 7^ and (m) (cG*2 + ^£3) = 0. In the following subsections we consider all the three 
cases one by one and will show that the above cases (i), (ii) and (in) correspond to three, 
two and one parameter Lie point symmetries, respectively. 



A. Three parameter symmetry (Case (i) A2 = 4, c 7^ and {C.G2 + AG?) 7^ 0) 

Solving Eqs. (1581) and (|60|) for A2 = 4 with c ^ we get A3 = and the explicit forms 
for the functions a, c and d turn out to be 

a = ai + ^-(^t 2 + d 2 t), c =^-(^t + d 2 ), d=(^t + d 2 ), (63) 

where a±, d\ and d 2 are three arbitrary (symmetry) parameters, which lead to a three pa- 
rameter Lie point symmetry group. Substituting A 2 = 4, A 3 = and the actual forms of G 2 
and G3 in Eq. (ISTj) and integrating, we obtain 

-3 



g(x) = gie-f fd *(\i + J e^ dx dx) , (64) 
where g\ is an integration constant. For this form of g Eq. §2§ reads 

x + f(x) x 2 + gi e~ / fdx (Xj. + J e* fdx dx^ ~ 3 = 0. (65) 
The infinitesimal symmetries associated with Eq. f )65|) are then 

e = ai + l(^ 2 + rf 2 t), V= ± e -Sf**(x 1 + J e^ dx dx){^t + d 2 ). (66) 

The corresponding infinitesimal generators can be written as 

g t 2 g te~S fdx (\ l + j e$ fdx dx^j g 

Xl = dt' X2 = 2\~ 1 dt + 2A^ 0? 

2t d e-Sf dx (\ 1 + je$f dx dx) g 

X3 ~ T l di + A^ d~x (67) 

This represents a three parameter symmetry case. The corresponding Lie algebra is 

[A 1; X 2 ] = ^A 3 , [A 1; A 3 ] = |-A 1; [A 2 , A 3 ] = -|-A 2 . (68) 
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Example: 



As a specific example of three parameter Lie point symmetries we consider the nonlinear 
ODE, with f(x) = -§ and g 1 — X in Eq. (jHSD, 

x--x 2 + - -^3=0. 69 

x (Ai x — iy 



The infinitesimal symmetries for this equation can be written from Eq. (1661) as 

t = a l + ^e + d 2 t), V=(^^-)(jt + d 2 ). (70) 
The associated infinitesimal generators read 

Xi = l' A ' 2 = ^-| + (^r 1 ) 4. * = r I + (^4^) < 71 > 

at 2Ai ot V 2 Ai J ox \\ at \ X\ J ox 



The general solution of Eq. (1691) is given by 22] 



_ Ai(mAi T >/mAf (m - jjAfj + m 3 (t - t ) 2 ) 
where I±, to are constants of integration and m = I±X 2 — A. 

B. Two parameter symmetry (Case (ii) A2 7^ 1, 4, c = and (cG2 + riG^) 7^ 0) 

Again solving Eqs. (J58l) and ( l60l) consistently for the choice A2 7^ 4, c = with (cG 2 + 
dGs) 7^ 0, we get A3 = and the explicit forms for the functions a, c and d as 

2a 2 2a 2 Ai 

a = ai + a 2 t, c = — , a = — - — , (73) 
A 2 A 2 

where ai and a 2 two arbitrary constants which corresponds to two parameter Lie point 
symmetry group. Substituting A3 = and the actual forms of G 2 and G3 in Eq. ( 16TT) and 
integrating, we obtain 

g(x) = g 1 e-f fdx (x 1 + j e^ dx dx) \ A 2 + 4, (74) 

so that Eq. (jSJ) takes the form 



x 



+ f(x)x 2 + g 1 e-f fdx (x 1 + J ef fdx dx^ * = 0, A 2 ^ 4. (75) 
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Note that in Eq. (|74p . for the choice A2 = 1 we have g(x) = g±e~ ■> f dx , where g\ is an 
arbitrary constant. This is exactly the same case ( 1381) with g\ = (#2 is an arbitrary 
constant which can be considered as g\ for analogy) admitting eight Lie point symmetries. 
So we exclude the A 2 = 1 case in Eq. (1751) from the present two parameter symmetry case. 
The infinitesimal symmetries associated with Eq. (175]) are 

Z = a 1 + a 2 t, r] = ^-e-S fdx (\ 1 + j e^ dx dx). (76) 

The corresponding infinitesimal generators become 

*l + ^- /J *(^ + / ef "« b ) I- < 77 » 

satisfying the commutation relation 

[X 1 ,X 2 ]=X 1 . (78) 

As an example for the two parameter symmetry group case, we consider the nonlinear 
ODE (for A 2 = 3) 

2 „ Ax 4 

x--x 2 + - — =0. 79 

x (Aix — iy 

The infinitesimal symmetries of Eq. ( 1791) can be written from Eq. (1761) as 

2(32 9 

£ = ai + a 2 £, r]=— (Aix — x). (80) 
Then infinitesimal generators can be written as 

The solution of (179"j) can be written in an implicit form 

A A^ 2 xln^ — Yj2~ ( m — ^! x ( m + ^) — ^/ w(Aix — l)(/iA^x — 1) j j 
A A-^ x 

-Jx^m^x - l)(hXjx -m)± m 3/2 (t - t Q )x = 0, (82) 
where I\ and to are integration constants and m = I\X\ — 2 A. 
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C. One parameter symmetry (cG2 + dG% = 0) 



Next we consider the case CG2 + dG^ = 0. As G2 7^ and G3 7^ we have necessarily 
c = and d = 0. With this set of restrictions in f[2"2"j) and f[2"31 we get d = 0, <? 7^ 0, which 
gives a = ai, where a\ is an arbitrary constant. Obviously this leads to an one parameter 
symmetry group. It corresponds to time translation generator X = irrespective of the 
form of / and g. Hence all other forms of / and g which do not belong to eight, three, and 
two parameter symmetry groups belong to the family of one parameter symmetry group. 

Now, multiplying with an integrating factor xe 2 ^^ dx , Eq. (T5]) can be written as 

£ e 2 //(*)«&: (x + f{x)i?) = -g(x)xe 2 f fix)dx (83) 
After an integration the above equation can be brought to the form 

rf e 2fm<h> + 2 J gfr) e 2jf(x)dx dx = (g4) 

where I\ is an integration constant. A second integration leads to the quadrature 

±(t-t )= [ dX (85) 

J J e -2ff(x)dx( Il _2fg( X ) e 2im*'dx) 

where to is the second integration constant. Depending on the form of f(x) and g(x), one 
may or may not be able to carry out the integration on the right hand side of (j85p explicitly. 
We now consider two specific examples belonging to this class. 



Example 1: Mathews- Lakshmanan Oscillator \22J 
The equation of motion is 

x Xx x 2 + U ° X = (86) 

2 

with f(x) = — 14 X A x r 2 and g(x) = j-^p- It admits only the translational symmetry as point 
symmetry. Using these forms of / and g in Eq. ( |85l) . we can write down the general solution 

as 



22] 



x(t) = Asm(nt + 6), n = - U ° (87) 
W V ; VI + \A 2 V ; 



Eq. f )86|) has a Hamiltonian 



2 ™2 



H = 5^ 1 + ^) + 5(iTIPy (88) 
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where the canonical conjugate momentum p = 1+ x Xx 2 ■ Note that when A is negative, | x |< 
The corresponding phase space structure for Eq. ( 1861) is given in Fig. 3. 




1.0 -0.5 0.0 0.5 1.0 



FIG. 3: (Color online) Phase portrait of equation (j86|) for different values of A corresponding to 
the Hamiltonian (I88p . (a) A > and (b) A < with the restriction | x |< 



Example 2: Particle in a rotating parabolic well ]2a] 



Similarly a particle moving in a rotating parabolic potential well having the equation of 
motion 



x + 



A x 



is a well known model 



22 



26 



Ax 2 



if + 



Uq X 



l + Xx 2 



0, A > 



271 ]. The underlying Hamiltonian is found to be 



1 



2 (1 + \x*) + 2 Uj2 ° x2 ' 



(89) 



(90) 



and the canonical momentum 



p = x (1 + A x 2 



(91) 



Eq. ( 1891) admits only the translational symmetry as the admissible Lie point symmetry. 
Using Eq. (|84p one can write the first integral for Eq. fl89l) as 



(1 + Ax 2 ) x 2 +w 2 i 2 = h. 



(92) 
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The above equation can be integrated to get the solution in terms of complicated elliptic 
integral of second kind. 

Considering the general equation corresponding to one parameter Lie point symmetry, as 
noted above, one can deduce the first (energy) integral from the time translational symmetry 
which in turn can be integrated to yield the general solution. One can also rewrite the first 
integral of the above such examples as a linear first-order equation in some new variables, 
which can be integrated to give the second integral explicitly. Now, the question arises why 



the above mentioned examples are integrab 
symmetry. Recently, Bhuvaneswari et al. 



e explicitly when they admit only one Lie point 



29| proved that Eq. ( 155]) exhibits the so called 



nonlocal A-symmetries as the possible reason for the integrability of such systems. Hence 
one can explain the integrability of such systems by associating more general symmetries to 
them. So there exists the question of identifying the systems belonging to the family of one 
parameter Lie point symmetry that are explicitly integrable. This problem has not been 
pursued in this study and will be discussed separately. 

VII. LIE SYMMETRIES OF EQ. © WITH f(x) = OR g(x) = 
If we consider the special case of Eq. ([2]) with f(x) = 0, then we have 

x + g(x) = 0. (93) 



The above equation has been studied in detail by Pandey et al. [12j and it has been shown 
that it admits two and three parameter symmetry generators for certain forms of g(x), while 
for other choices it admits only one parameter symmetry group. 
Next, for the case g(x) = 0, Eq. (T5]) takes the form 

x + f(x)x 2 = 0. (94) 

To get the symmetry generators of Eq. ( 1941) we substitute g as zero in Eqs. ( |T9i) and ( 1201) and 
Gix, Gix and G% x in terms of G\, G2 and G3 and then equate the coefficients of G±, Gi and 
G3 to zero. Doing so we arrive at the following determining equations for the infinitesimals, 
that is 

6 = 0, c = 0, d = 0, a -26 = 0. (95) 
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Solving these equations one can get the forms of a, b, c and d which on substituting in Eq. 
(TIT]) and ( !I8|) give the symmetry generators as 

£ = ai + a 2 t + cit 2 + (bit + 6 2 ) 

?y = 2 bxG x + (c x t + c 2 )G 2 + (d x t + d 2 )G 3 , (96) 

where 3, G\, G 2 and G3 are already defined in Sec. II. Here ai, a 2 , 61, 62, ci, c 2 , a?i and d 2 are 
the eight symmetry generators. 

The corresponding infinitesimal generators can be then written as 

d d d d d 

2 d 

X5 — t -57 + £^2— , ^6 = C 2 — , X 7 = tG 3 — , X 8 = G 3 — , (97) 
at ox ox ox ox 

satisfying the s/(3, R) algebra. Obviously Eq. ( |94|) can be linearized under the transforma- 
tion X = h(x) to the free particle oscillator equation 

1 = 0, (98) 

where h(x) = hi J e^^ dx dx + h 2l hi and h 2 being arbitrary constants. 



VIII. EQUIVALENCE TRANSFORMATIONS 



Finally we consider equivalence transformations 3lj associated with ([2]). Let us consider 
a set of smooth, locally one-to-one transformations T : (t,x,f,g) — > (T, X, fi,gi) of the 
space R 4 that act by the formulae 



T = G(t,x), X = F(t,x), fi = H(t,x, /), gi = L(t,x,g). 



(99) 



An equivalent transformation of Eq. (J2J) is an invertible transformation that converts Eq. 
(J2]) to an equation of the same form 31 1 



X = -fi{X)X 2 -gi{X). 



(100) 



In this case Eqs. ([2]) and (11001) and the functions {f(x),g(x)} and {fi(X),gi(X)} are equiv- 
alent. 
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Substituting the transformation (199]) into Eq. (11001) we get 

fx(F t + xF x ) 2 (G t + xG x ) + 9l (G t + xG x f = (F t + xF x )[(G tt + 2xG tx 

+ x 2 G xx - G x (fx 2 + g)] - (G t + ±G x )[{F tt + 2xF tx + x 2 F xx - F x (fx 2 + g)], (101) 

where the subscripts denote partial derivative with respect to that variable. Equating the 
coefficients of different powers of x n , n — 0, 1, 2, 3, we get 

f\F 2 G x + giG x = F X G XX — G X F XX , (102) 

h(2F x F t G x + F 2 G t ) + 3 9l GlG t = Jf + F t G xx - G t F xx + 2F x G tx - 2G x F tX} (103) 

h(2F x F t G t + F 2 G X ) + 3 9l G 2 G x = F x G tt - G x F tt + 2F t G tx - 2G t F tx , (104) 

hF 2 G t + gi G] = G tt F t - G t F tt + gj, (105) 

where J = G t F x — F t G x . Solving Eqs. f 1 1 2 j) and (11031) we get the forms of the functions fi 
and g\ as 

2GtG x F xx — 3F x GtG xx + fG x J + F t G x G xx + 2F x G x Gt x — 2G 2 F tx / 1r .„x 

/l = 2F-G7J ' (1 ° 6) 

F x G xx (G t F x + F t G x ) + F x G x (2(G x F tx — F x G tx ) — fJ) — 2G 2 x F t F xx , in _s 
9i- • ( 107 ) 

Substituting these values of fi and g\ in Eqs. f ll04j) and H105[) and then equating both the 
equations thus obtained, we get a general condition on the forms of the functions F and G 

as 

{Gt — 1)[(«/ + 2f t G x ){2G 2 x F tx + JG XX ) — G x (2G x (F t G t F xx + F x G x F tt ) + F x G t (2F x G tx 
+/^))] + (-ft^x — F X )[2G 2 F X G XX + 2G 2 X F X G U — 2F x G x G t G tx — fJG x G t ] + J[2F x G x G tx 
+f(JG x + 2G^F t ) - 2^F ffi C^] = 0. (108) 

The above equation is difficult to solve for the general forms of -F and G. Hence one can 
consider specific forms of F and G for simplicity. One of the possible solutions for the set 
of Eqs. flil)2>(fI05]) is of the form 

F = ax + /?, G = jt + S, (109) 

where a, (3, 7 and 5 are arbitrary constants. Substituting the above forms of F and G in 
Eqs. (1102p - (ll05j) . we get the forms for the functions fi and g\ as 

h = ~, 9i = -5- ■ 110 
22 



Then the possible equivalence transformation is 



X = ax + P, G = -ft + 5, (111) 

fi = ~, 9i = -r ■ (112) 

a 7 Z 

Making use of the above, we can write down the equivalence transformation for some of the 
examples discussed earlier. For example, for Eq. f )69|) the transformation is given as 

X = ax + (3, G = jt + S, (113) 

2_ A(X-/3) 5 

h X - ' 91 a 7 2 (A!(X - /3) - a) 3 ' lii4j 

whereas for Eq. (1791) the transformation is 

X = ax + (3, G = jt + 5, (115) 
f - 2 A(X - (5 f 

h ~ "x^fl' 91 - <ry*(Xi(x - p) - <*)* ■ (116) 

One can analyse Eq. fl 1 8 [) further to get more general equivalence transformations. We 
have not pursued this problem further in the present work. 



IX. CONCLUSION 



In this paper we have investigated systematically the Lie point symmetry groups asso- 
ciated with the quadratic Lienard Eq. (j2j). Even though the integrability properties of 
the general form of Eq. (j2]) has been discussed to some extent in the literature, we have 
systematically identified and classified all those equations which admit one, two, three and 
eight parameter symmetry groups. We have found the general form of (J2]) that belongs to 
the linearizable case as the one given by Eq. (|25|) admitting the maximal (eight parameter) 
symmetry group, whereas for the integrable but not linearizable cases the general forms for 
three and two parameter symmetry groups are represented by Eqs. (165j) and ( 1751) . respec- 
tively. We have also deduced the interesting result that the condition for isochronicity of 
Eq. 02]) is the same as that of the linearizability condition. Our analysis clearly confirm the 
powerful nature of analysis of nonlinear ODEs based on symmetry properties. 



23 



X. ACKNOWLEDGMENTS 



AKT and SNP are grateful to the Centre for Nonlinear Dynamics, Bharathidasan Univer- 
sity, Tiruchirappalli, for warm hospitality. The work of SNP forms part of a Department of 
Science and Technology, Government of India, sponsored research project. The work of MS 
forms part of a research project sponsored by UGC The work forms part of a Department of 
Science and Technology, Govt, of India IRHPA project and a Ramanna Fellowship project 
of ML. He also acknowledges the financial support provided through a DAE Raja Ramanna 
Fellowship. 



[1] N. H. Ibragimov, Elementary Lie Group Analysis and Ordinary Differential Equations (Wiley, 
New York, 1999). 

[2] G. Baumann, Symmetry Analysis of Differential Equations with Mathematica (Springer- 
Verlag, Berlin, 1998). 

[3] G. W. Bluman and S. C. Anco, Symmetry and Integration Methods for Differential Equations 

( Springer- Ver lag, New York, 2002). 
[4] P. Hydon, Symmetry Methods for Differential Equations: A Beginner's Guide (Cambridge 

University Press, Cambridge, 2000). 
[5] J. M. Hill, Differential Equations and Group Methods for Scientists and Engineers (CRC 

Press, Boca Raton, 1999). 
[6] B. J. Cantwell, Introduction to Symmetry Analysis (Cambridge University Press, Cambridge, 

2000). 

[7] N. Euler and W. H. Steeb, Continuous Symmetries, Lie Algebras and Differential Equations 

(B.I. Wissenschafts verlag, Mannheim, 1993). 
[8] P. J. Olver, Applications of Lie Groups to Differential Equations (Springer- Ver lag, New York, 

1986). 

[9] H. Stephani, Differential Equations: Their Solutions Using Symmetries (Cambridge University 
Press, Cambridge, 1989). 
[10] M. Lakshmanan and S. Rajasekar, Nonlinear Dynamics: Integrability, Chaos and Patterns 
(Springer- Ver lag, Berlin, 2003). 



24 



[11] F. M. Mahomed, Math. Meth. Appl. Sci. 30, 1995 (2007). 

[12] S. N. Pandey, P. S. Bindu, M. Senthilvelan and M. Lakshmanan, J. Math. Phys. 50, 082702 

(2009); 50, 102701 (2009). 
[13] F. M. Mahomed and P. G. L. Leach, J. Math. Phys. 30, 2770 (1989). 
[14] P. G. L. Leach, J. Math. Phys. 22, 679 (1981). 

[15] P. G. L. Leach and V. M. Gorringe, J. Phys. A: Math. Gen. 23, 2765 (1990). 

[16] K. S. Govinder and P. G. L. Leach, J. Phys. A: Math. Gen. 27, 4153 (1994). 

[17] O. Gat, J. Math. Phys. 33, 2966 (1992). 

[18] A. R. Chouikha, J. Math. Anal. Appl. 331, 358 (2007). 

[19] M. Sabatini, J. Diff. Eq. 196, 151 (2004). 

[20] I. Boussaada, A. R. Chouikha and J. M. Strelcyn, Bull. Sci. Math. 135, 89 (2011). 
[21] M. Bardet, I. Boussaada, A. R. Chouikha and J. M. Strelcyn, Bull. Sci. Math. 135, 230 (2011). 
[22] P. M. Mathews and M. Lakshmanan, Quart. Appl. Math. 32, 215 (1974). 
[23] J. F. Carinena, M. F. Ranada, M. Santander and M. Senthilvelan, Nonlinearity 17, 1941 
(2004). 

[24] J. F. Carinena, M. F. Ranada and M. Santander, Rep. Math. Phys. 54, 285 (2004). 

[25] H Goldstein, Classical Mechanics (Narosa Publications, New Delhi, 2001). 

[26] A. Venkatesan and M. Lakshmanan, Phys. Rev. E 55, 5134 (1997). 

[27] A. H. Nayfeh and D. T. Mook, Nonlinear Oscillations (Wiley, New York, 1979). 

[28] C. Muriel and J. L. Romero, IMA J. Appl. Math. 66, 111 (2001); J. Lie Theory 13, 167 

(2003); J. Phys. A: Math. Theor. 42, 365207 (2009). 
[29] A. Bhuvaneswari, V. K. Chandrasekar, M. Senthilvelan and M. Lakshmanan, J Math. Phys. 

53, 073504 (2012). 

[30] F. Calogero, Isochronous Systems (Oxford University Press, USA, 2008). 
[31] L. V. Ovsyannikov, J. Appl. Mech. Tech. Phys. 45, 153 (2004). 



25 



